Astronomy & Astrophysics manuscript no. 7653 


© ESO 2008 


February 5, 2008 





o 
o 

(N 



On magnetic instabilities and dynamo action in stellar radiation 

zones 

J.-P. Zahn 1 , A. S. Brun 21 and S. Mathis 3 ' 2 ' 1 



i— H ' ' LUTH, Observatoire de Paris, CNRS, Universite Paris-Diderot; Place Jules Janssen, F-92195 Meudon, France 
^ , e-mail: jean-paul . zahn@obspm. fr 
^ ■ 2 DSM/DAPNIA/Service d'Astrophysique, CEA Saclay, F-91191 Gif-sur-Yvette, France; AIM, UMR 7158, CEA - CNRS - Universite Paris 7 
Observatoire de Geneve, 51 chemin des maillettes, CH-1290 Sauvemy, Switzerland 
e-mail: sacha.brun@cea.fr; stephane.mathis@cea.fr 



43 ■ 

Received 17 April 2007; accepted 18 July 2007 
8 . ABSTRACT 

r$ Context. We examine the MHD instabilities arising in the radiation zone of a differentially rotating star, in which a poloidal field of fossil origin 
is sheared into a toroidal field. 

Aims. We focus on the non-axisymmetric instability that affects the toroidal magnetic field in a rotating star, which was first studied by Pitts 
and Tayler in the non-dissipative limit. If such an instability were able to mix the stellar material, it could have an impact on the evolution of 
the star. According to Spruit, it could also drive a dynamo. 

Methods. We compare the numerical solutions built with the 3-dimensional ASH code with the predictions drawn from an analytical study of 
the Pitts & Tayler instability. 

£C) ' Results. The Pitts & Tayler instability is manifestly present in our simulations, with its conspicuous m = 1 dependence in azimuth. But its 
analytic treatment used so far is too simplified to be applied to the real stellar situation. Although the instability generated field reaches an 
energy comparable to that of the mean poloidal field, that field seems unaffected by the instability: it undergoes Ohmic decline, and is neither 
eroded nor regenerated by the instability. The toroidal field is produced by shearing the poloidal field and it draws its energy from the differential 
, rotation. The small scale motions behave as Alfven waves; they cause negligible eddy-diffusivity and contribute little to the net transport of 
• • angular momentum. 

, Conclusions. In our simulations we observe no sign of dynamo action, of either mean field or fluctuation type, up to a magnetic Reynolds 
, number of 10 5 . However the Pitts & Tayler instability is sustained as long as the differential rotation acting on the poloidal field is able to 
■ generate a toroidal field of sufficient strength. But in the Sun such a poloidal field of fossil origin is ruled out by the nearly uniform rotation of 
' the deep interior. 
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1. Introduction 

In recent years there has been a surge of interest for stellar mag- 
netism, due mainly to the discovery of magnetic fields in an in- 
creasing number of stars, and to their mapping through Zeeman 
imaging (cf. Donati et al. 2006). Theory benefits enormously 
from these new constraints, and quite naturally the main fo- 
cus is on the generation of magnetic fields in turbulent convec- 
tion zones, which can now been studied through high resolu- 
tion numerical simulations (Brun et al. 2004, 2005; Dobler et 
al. 2006). But attention has been drawn also to the instabilities 
that may affect the magnetic field in stably stratified radiation 
zones. Spruit (1999) reviewed various types of instabilities that 
are likely to intervene in a magnetized radiation zone, and he 
concluded that the strongest among them were those which had 
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been described by Tayler and his collaborators. Indeed, Markey 
and Tayler (1973) have shown that a purely poloidal field would 
be unstable to non-axisymmetric perturbations, and so would 
also a toroidal field (Tayler 1973; Wright 1973; Goossens et 
al. 1981). Later Pitts and Tayler (1985) proved that even in the 
presence of rotation a toroidal field would be unstable to such 
perturbations. 

Spruit (1999, 2002) analyzed the latter instability in more 
detail, including Ohmic dissipation and radiative damping, and 
starting from the dispersion relation established by Acheson 
(1978); he suggested that it could regenerate the toroidal 
field, and thus drive a genuine dynamo. Applying essentially 
Spruit's prescriptions, Maeder and Meynet (2003, 2004, 2005), 
Eggenberger et al. (2005) and Heger et al. (2005) introduced 
this dynamo, and the turbulent transport believed to be asso- 
ciated with it, in their stellar evolution calculations. It proved 
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extremely efficient in establishing quasi uniform rotation. In 
the case of massive stars, it would increase the loss of angular 
momentum through the stellar wind, and lead to slower rotat- 
ing stellar cores; this would yield pulsar rotation rates that are 
in better agreement with the observations. 

However, guided by other observational evidence, 
Denissenkov and Pinsonneault (2006) questioned the existence 
of such a powerful process of angular momentum transport; 
they noticed some inconsistencies in Spruit's heuristic ar- 
gumentation, which they tried to correct. Unfortunately the 
choice they made was in conflict with the rigorous analytical 
treatment of the instability; this was pointed out by Spruit 
(2006), and that part of their discussion was deleted in the final 
version (Denissenkov & Pinsonneault 2007). In the mean- 
while, Braithwaite (2006) conducted a numerical experiment 
which he claimed to validate Spruit's dynamo scenario. 

Our interest in those instabilities was aroused when we ob- 
served them while verifying if a fossil field was able to pre- 
vent the spread of the solar tachocline, as had been proposed 
by Gough and Mclntyre (1998): the results of these three- 
dimensional simulations are reported in Brun and Zahn (2006). 
We then wanted to check whether these calculations agreed 
with Spruit's analytical predictions, and this incited us to re- 
examine his original derivation. 

We begin by extending Acheson's dispersion relation to the 
case where both an entropy gradient and a composition gra- 
dient are present, and derive the actual solutions of that equa- 
tion (rather than upper or lower bounds). We then discuss the 
heuristic arguments that were first used by Spruit, and show 
how they can be misleading. Next we compare the analytical 
solutions with the numerical results obtained in Brun & Zahn 
(2006, hereafter referred to as BZ06), and finally we examine 
whether the dynamo that was suggested by Spruit does actually 
operate. 

2. The non-dissipative case 

Let us first recall the major results concerning the instability 
that affects a toroidal field in the radiation zone of a rotating 
star, which were derived by Pitts & Tayler (1985) in the non- 
dissipative limit. We consider this field as given, and ignore 
what causes it, namely the mean poloidal field and the differen- 
tial rotation that is acting on it; we further assume that it varies 
as By cc m p with the distance from the rotation axis (p > 1). To 
ease the task for the reader, we put the technical developments 
in Appendix. 

As was done in previous works, we perform a local analysis 
and examine the stability of an imposed axisymmetric toroidal 
field. We submit it to wave-like perturbations of the form 



exp i(lm + imp + nz — crt), 



(1) 



in a cylindrical reference frame (zzr, <p, z) centered on the rota- 
tion axis of a star that is rotating uniformly with angular ve- 
locity Q; cr = <t r + icrj is the complex frequency. Such per- 
turbations obey a dispersion relation that was first derived by 
Acheson (1978, Eq. 3.20); Spruit (1999) made it more tractable 
by applying it only to the vicinity of the rotation axis, by assum- 
ing that (//«) <k 1 and that the Alfven velocity is negligible 



compared to the sound speed. In Appendix, we establish that 
equation in a slightly more general case, allowing the stratifi- 
cation to depend both on entropy and on chemical composition. 
For m + and ignoring all forms of energy dissipation (radia- 
tive, Ohmic, viscous) the dispersion relation takes the form (cf. 
Spruit 1999, Eq.A9) 



(cr 2 - co 2 A ) (cr 2 - u\[l+ A]) + (p- l)^f (imQcr + 2cr 2 ) 



2Q.O- + 



p+l 



20cr + 2- 



(2) 



We have introduced here the stratification parameter 



I 2 N 

A = — — 

n L a>2 



(3) 



which measures the relative strength of the two restoring 
forces: buoyancy and magnetic tension; N is the buoyancy fre- 
quency (cf.[8]) and ioa the Alfven frequency associated with the 
toroidal field B x : 



(¥) 



m 2 B 2 
m 2 Anp 



(4) 



This dispersion equation admits two pairs of wave modes, 
that are well separated in frequency when u>a «c £2, which 
is the case we shall consider here of moderate field strength 
(By <K 100 kG below the solar convection zone). The fast 
modes may be isolated by setting cr 2 » a>\, which yields to 
first approximation cr 2 - (l/n) 2 N 2 + 4Q. 2 : we recognize here 
the signature of the gravito-inertial waves. The slow modes 
are found in the domain cr 2 <K <y 2 <K Q 2 , where the dis- 
persion relation reduces to the following quadratic equation in 
( = 2mQ.cr I u 2 A 



?+<\(- [m 2 (l + A) - 2(p + 1)] = 0. 



(5) 



Its discriminant is A' = m 2 (l + A) + 2(1 - p); when A' < 0, it 
has a complex root with a positive imaginary part 07 > and 
hence this mode would be unstable for 



p > I + —(1 + A). 



(6) 



As Pitts and Tayler (1985) pointed out, all modes are stable for 
p = 1 . And one expects actually the toroidal field to vary lin- 
early with m in the vicinity of the rotation axis, if it is generated 
by shearing a poloidal field through a depth dependent rotation. 
Hence we shall assume from here on that p — 1 . We are then 
left with two oscillatory modes, whose frequencies are given 
by 



-2±m Vl + A or cr = — 

£2 



1 VTTa 
— ± 

m 2 



(7) 



Note that for A < 3 there are two negative roots, and that for 
A > 3 the roots have opposite sign. 
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3. Overstable modes 

Since all modes are stable for p = 1 in the non-dissipative 
limit, instability can only occur in form of a diffusive instabil- 
ity, with the complex frequency presenting a positive imaginary 
part. The dissipation is provided here by radiative and Ohmic 
diffusion, with respective diffusivities k and rj. In addition to 
the two slow oscillatory modes discussed above, whose ampli- 
tude may then grow exponentially in time, there is also a direct 
mode, whose frequency cr R — > in the non-dissipative limit, 
that may become unstable. 

We consider here only the case where the Roberts number 
e = tj/k may be treated as a small quantity, which is the case in 
stellar interiors (s « 10 4 below the solar convection zone). In 
Appendix, we derive Acheson's dispersion relation in the more 
general case where both a gradient of entropy and of chemical 
composition are present. Then the stratification is characterized 
by two buoyancy frequencies: 



N 2 = N 2 + N 2 = — 
' >* H f 



(V ad - V) + 



g (dlnfi 



Hp \dlnP 



(8) 



with the usual notations, and assuming the perfect gas law for 
simplicity. The corresponding stratification parameters are then 



I 2 N 2 

A, = —-j andA^ — 

n l oj- n l oj a 



I 2 N 2 



(9) 



it is convenient to introduce also the sum A* = sA, + Au- 
la Appendix (Eq. IA. 15b we show that the scaled com- 
plex frequency £ = 2mQ. cr/co^ obeys the following third-order 
equation: 



sA,f + df+A^-A^S 2 

- is <r(<r 2 + 4<r + 3- a* -s 2 ) = o, 



(10) 



where S = 2Q.rjn 2 /u> A represents the suitably scaled Ohmic 
diffusivity. Note that A, intervenes only through the product 
sA t , which illustrates how the restoring force due to the en- 
tropy gradient is weakened through radiative damping. For 
given meridional wavenumbers n and /, or equivalently for a 
set [A,, A^, 5], Eq.[l0]has three complex solutions; if its imag- 
inary part is positive, the mode is overstable. Its growth- 
rate ctj = (a> A /2mQ.) d then peaks at azimuthal wavenumber 
\m\ — 1, a distinctive property of this instability. From here on 
we shall consider only that case \m\ = 1, and refer to this insta- 
bility as the Pitts & Tayler instability, to distinguish it from that 
which affects a poloidal field (Markey & Tayler 1973). 

We observe that if S , sA, and/or A^, are of order unity, this 
will also be the case for £, which implies that both the fre- 
quency cr R and the growth-rate 07 scale as (Ji/Cl, like in the 
non-dissipative case. The maximum growth-rate and the corre- 
sponding frequency are given in TableQ]for a sample of values 
of eAt/Ap, we see that it is the prograde mode that is unstable 
(or rather overstable), except for A^ = 0. 

For given ratio A^A^, the instability domain can be de- 
lineated in the [A*, S] plane by locating the points where the 
growth-rate vanishes and thus where the oscillation frequency 
is purely real: cr - cr R . Referring back to ( TTOb . this occurs when 



the following equations are both satisfied: 
eA,x 2 + (2x 2 + Ax-A ll )S 2 = 0, 



x 2 +Ax + 3-A* -S 2 



0. 



(11) 



where x = 2mQ. cr R /a> A is the scaled frequency. The solution of 
this system determines the minimum strength of the magnetic 
field that is required for instability, as we shall see next. 

Let us first examine the two limit cases (A, = and A^ = 0) 
that were considered by Spruit. 



3.1. Stratification due only to the composition gradient 

For A, = 0, system (fTTT i reduces to 

2x 2 +Ax-A ll =0, 
x 2 +Ax + 3-A ll -S 2 = 0. (12) 
It is easily solved: 

x = -1 + Jl +A M /2 

S 2 = 2 A /l+A /J /2+l-A jU /2; (13) 



the function S (A^) decreases from V3 at A^ — to at A max = 
6 + 4 a/3 = 12.928, as illustrated in Fig.Q] 

We now use this result to determine the instability thresh- 
old. Like Spruit (1999), we eliminate n 2 by combining the def- 



inition of A^: 

— JL- a 

2 2 — t 1 
n u 2 A 

with that of S : 
2 -^n 2 = S(A,); 



(14) 



(15) 



this gives us the minimum strength of the toroidal field required 
to trigger the instability: 



(16) 



[A M 5(A /J )]^=2Q ?7 /^. 

That condition is optimized when [A^SiA^)] is maximum, at 
[A5] max = 9.829, which is reached for A^, = 8.788. Thus our 
result confirms that of Spruit (1999, Eq. A23), although they 
differ somewhat in their numerical coefficients (A max , for in- 
stance, is underestimated by Spruit, whose value is 3 instead of 
12.928). 

3.2. Stratification due only to the entropy gradient 

For A A , = system (fTTT i has no simple analytical solution, but it 
can be easily shown that S (A,) vanishes at sA t = and 3, as il- 
lustrated in Fig.Q] Thus the maximum value of the stratification 
parameter that allows for instability is A, = 3/e. 

Proceeding as before, we find that the instability threshold 
is given by 



[eA t S(eA,)]oA 



2D,ET]l 2 N 2 ; 



(17) 



it is optimized at [sA,S(sA,)] = 0.728 for sA, = 2.070. Here 
also we retrieve Spruit's result, but with the exact proportion- 
ality constant. 
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Fig. 1. Limits of the instability domain in the (A*,S) plane, 
for various values of A t /A^. The stratification parameters are 
defined as A, = (l/n) 2 (N,/oj A ) 2 and A M = (7/n) 2 (Ay« A ) 2 , 
whereas A* = eAt+A^. s = tj/k designates the Roberts number, 
and S = 2rj£ln 2 l<x> 2 A defines the limits of the instability domain 
in vertical wavenumber n, for given A*. 

Table 1. Maximum scaled growth-rate £]™ x = 2£2o7/w 2 and 
corresponding frequency £™ x = 2Q.ct r I<jl> 2 , A* and S , at given 



and for 


m\ = 1. 










A* 


5 (A*) 


wnax 
<>« 


wliax 





3.306 


0.3463 


0.3868 


0.1722 


0.01 


3.203 


0.3593 


0.3736 


0.1679 


0.1 


2.957 


0.3774 


0.2820 


0.1397 


1 


1.817 


0.3750 


0.1160 


0.06945 


10 


1.908 


0.2599 


0.02211 


0.01736 


oo 








-0.4858 


0.4227 



3.3. Stratification due to both entropy and composition 
gradients 

In the general case both causes of stratification are present, 
and S is a function of A, and A^; the instability domain is dis- 
played in Fig. [TJfor a set of fixed A,/A M . All solutions, except 



for Aft = 0, have in common the point A* = 0, S 



V3, and 

all solutions, except for A, = share the point A* = 3, S = 0, 
as predicted by (fTTI) . Note that they are double-valued in the 
general case, which means that the instability interval in «, the 
vertical wavenumber, has a lower boundary that differs from 
zero. Only in the limit cases, A, — 0, A^ = 0, does the instabil- 
ity domain extend to S = 0, between A* = and A* max . 

To determine the instability threshold for the toroidal field, 
we eliminate as before n 2 between 



I 2 



N 2 
f 



= A* 



(18) 



Table 2. Stratification parameters and optimal conditions for 
the onset of instability. A* max is the highest value of A* = 
eAr+Aft that allows for instability; [AS ] mdx the maximum of the 
product A* S (A*) for fixed sAt/A^, obtained at A* pt , from which 
one derives the instability threshold for the magnetic field (see 
text). 



eAt/A^ 



S(KJ [AS] n 



A* 






12.928 





9.829 


8.788 


0.01 


10.249 


0.602 


9.307 


8.138 


0.1 


6.970 


0.754 


6.960 


5.885 


1 


4.000 


0.550 


3.451 


3.171 


10 


3.125 


0.219 


2.186 


2.087 


oo 


3 





0.728 


2.070 



and 



2Q.n 7 

— n 

oj 2 a 



s (A*y, 



this yields 

[A* S(A*)]lo 4 a = 2Q,T]l 2 [sN? 



(19) 



(20) 



which validates Spruit's prescription of replacing the buoyancy 
frequency A^ 2 by [sN? + Af 2 ] (his Eq. 55)0 This expression is 
optimized for the values of [A* S (A*)] that are given in Table|2] 
together with the maximum value of the stratification param- 
eter A* nax which allows for instability, and that of the corre- 
sponding S (A* max ). 

For either s = or N 2 = 0, we retrieve the two special 
cases considered before. 



4. On the validity of heuristic arguments 

In his seminal paper, Spruit (1999) used a few heuristic argu- 
ments to describe the Pitts & Tayler instability in a simple, intu- 
itive way. Let us examine their validity in the light of the rigor- 
ous treatment given in the preceding section and in Appendix. 

With one such arguments, Spruit seeks a lower limit for the 
vertical wavenumber n. We quote him: "For displacements of 
amplitude the work done per unit mass against the stable 
stratification is j£ 2 (l/n) 2 N 2 . The energy gained from the field 
configuration is jcj a ^ 2 . For instability, the field must be strong 
enough, such that oj 2 a > (l/n) 2 N 2 " Furthermore, since I 1/r, 
this provides him a lower limit for n (cf. his Eq. 44). 

In Spruit (2002), the argument is slightly different: £ is the 
unstable displacement, and jw 2 ^ 2 the kinetic energy that is re- 
leased in this displacement, "since the growth-rate of the insta- 
bility in the absence of constraints like rotation is <x ~ Oa". 

Another argument is used to set an upper limit to n: the 
growth-rate <x has to overcome the damping rate of Ohmic dif- 
fusion; hence 

cr>n 2 t], (21) 
and this is applied to the case Q » oja, where 



1 Through heuristic arguments, Maeder & Meynet (2004) obtain a 
similar expression, but where A', 2 is reduced by e/2 instead of e. 
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Denissenkov & Pinsonneault (2006) object quite under- 
standably to this ambivalent choice of cr, namely cr ~ u> A for 
the lower bound and cr ~ io 2 A /Q for the upper bound on n. They 
take the latter definition in both cases, which seems more co- 
herent - but disagrees with the rigorous treatment. 

Where is the problem? In fact buoyancy and magnetic field 
act both as restoring forces, as is illustrated by the energy equa- 
tion, which we establish in Appendix in the non-dissipative 
case ( 1A.221 i: 



U d 4\ + 1 -mW A f + 1 -N 2 £=0. 
2\dt 2 A 2 z 



(22) 



It shows how the energy of the perturbation is exchanged back 
and forth between kinetic and potential (magnetic and buoy- 
ancy) energy, as is typical in dynamical systems. In the ab- 
sence of dissipation, the total energy remains constant, while it 
increases when the diffusive instability occurs. Thus the buoy- 
ancy force and the Lorentz forcdl are not antagonistic, and it is 
not possible to draw an instability condition from their relative 
strengths. 

These shortcomings of the heuristic approach have been 
implicitly recognized in Spruit (2006), where he writes that 
"all physical effects contribute if all terms and factors [meaning 
A^, sA t , S and £ in our notation] are of the same order [i.e. of 
order unity, as can be seen in (fTUlil". To achieve this, the ratio 
{I In) 2 has thus to satisfy {I In) 2 ~ io 2 A IN 2 , a condition that is 
less stringent than Spruit's original inequality. 

The other heuristic condition ( |2TT i deserves also a couple of 
remarks. First, as we show in Appendix, in the limit of small 
n 2 the growth-rate is proportional to the damping rate: cr oc n 2 rj 
(cf. IA.16l l. and it is not clear what information can be extracted 
then about n 2 by equating these two rates. Second, in presence 
of radiative damping, why should this condition not be replaced 
by cr > n 2 K, since the thermal diffusivity k largely exceeds the 
Ohmic diffusivity? 

We conclude that heuristic arguments, though appealing, 
can be misleading when they are not supported by a rigorous 
treatment. 



5. Can the Pitts & Tayler instability sustain a 
dynamo? 

We now come to the suggestion made by Spruit (2002), that this 
instability could sustain a dynamo in stellar radiation zones, as 
thermal convection does in a convection zone. We find his idea 
quite interesting, but we argue that this dynamo cannot operate 
as he describes it. According to him, the instability-generated 
small scale field, which has zero average, is wound up by the 
differential rotation "into a new contribution to the azimuthal 
field. This again is unstable, thus closing the dynamo loop." 
But this shear induced azimuthal field has the same azimuthal 
wavenumber as the instability-generated field, i.e. m + and 



2 In fact, one should call it the Laplace force, contrary to what has 
become common practice. Pierre-Simon Laplace (1749-1827) was the 
first to give the expression of the force exerted by a magnetic field 
on an electric current; the force named after Hendrik Lorentz (1853- 
1928) is that experienced by a moving charged particle. 



Only possible dynamo loop: Vx < v x b > 



Axisym metric 
poloidal field 



B m (m -O) 



X. 



Q-effect 

Vft (m = 0) 



High m *0 
instability 



Axisymmetric 
toroidal field 



B T (m =0) 



X 



Taylor's 
instability 



Non-axisymmetric 
instability-generated field 



b (m * 0) 



Loop A 



Q-effectTO (m = 0) 
- leads to b(m^O) not to B(m=0) 



Loop B 

Fig. 2. How to close the dynamo loop involving the Pitts 
&Tayler instability. In dashed lines the loops proposed by 
Spruit (A) and Braithwaite (B). The only possible way to re- 
generate the mean toroidal and/or poloidal field is through the 
mean electromotive force < v x b > produced by the non- 
axisymmetric instability-generated field. But the dynamo must 
regenerate the poloidal field, in order to be fed from the dif- 
ferential rotation, and this leaves only the loop drawn in solid 
line. 

predominantly m = 1: it has no mean azimuthal component, 
and thus it cannot regenerate the mean toroidal field that is 
required to sustain the instability. For the same reason, the 
instability-generated field cannot regenerate the mean poloidal 
field, as was suggested by Braithwaite (2006). Therefore the 
Pitts & Tayler instability cannot be the cause of a dynamo, as 
it was described by Spruit and Braithwaite. 

To close the dynamo loop, as is well known in mean-field 
theory (Parker 1955; Moffatt 1978), one has to invoke the so- 
called a-effect, which involves the non-zero mean electromo- 
tive force < v x b > that is produced here by the Pitts & Tayler 
instability. This can be read in the azimuthally averaged induc- 
tion equation: 

dB 

— = e v {mB m ■ VQ] + V x < v x b > -V x (77V x B), (23) 

where v and b are the non-axi-symmetric parts of the veloc- 
ity and magnetic fields; the meridional advection has been ab- 
sorbed in the Lagrangian time derivative. The first term on 
the RHS describes how the poloidal field B m is wound up by 
the differential rotation to produce the mean toroidal field (the 
Q-effect), the second how the mean electromotive force may 
(re)generate both the poloidal and the toroidal fields (the a- 
effect), and the last term represents the Ohmic diffusion, even- 
tually enhanced by the turbulence (the j6-effect). The only pos- 
sible dynamo loop is depicted in Fig.|2](solid lines); those pro- 
posed by Spruit and Braithwaite are shown in dashed lines. 
Moreover, for the dynamo to operate, the mean electromotive 
force must overcome the Ohmic dissipation. 

Another type of dynamo has been observed in numeri- 
cal simulations, with externally forced turbulence: it is the 
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Parameter 


Symbol 


Sun 


Case A 


Case B 


thermal diffusivity 


K 


10 7 


810 12 


810 12 


magnetic diffusivity 


1 


10 3 


810 10 


810 9 


viscosity 


V 


30 


810 9 


8 10 9 


buoyancy frequency 


N, 


2.1 10" 3 


3 10- 4 


310- 4 


rotation frequency 


Q. 


3 10- 6 


3 10-" 


310- 6 



Table 3. Typical values of the relevant parameters in the upper 
radiation zone of the Sun, and values adopted for the numerical 
simulations (in cgs units). 



small-scale or fluctuation dynamo. There the magnetic field 
has no mean components, neither poloidal nor toroidal, but a 
quasi-stationary regime may be achieved when the magnetic 
Reynolds number Rm=V L/rj exceeds some critical value (V 
and L are the characteristic velocity and length-scale of the tur- 
bulence). But it remains to be checked whether such a dynamo 
can be sustained by an imposed shear, such as a differential 
rotation. 

6. 3-D numerical simulations of the MHD 
instabilities 

In Spruit's scenario, the toroidal field which undergoes the Pitts 
& Tayler instability is generated by winding up an existing 
poloidal field through a differential rotation Q(z) that results 
from slowing down the star by a stellar wind. It is similar to the 
model we took to examine the possibility of confining the solar 
tachocline by a fossil field (BZ06). However, in our case the 
differential rotation is imposed in latitude by the adjacent con- 
vection zone, and the depth dependence of Q near the polar axis 
is caused by thermal diffusion (cf. Spiegel & Zahn 1992). What 
distinguishes our model from that of Spruit, and several others 
(cf. Miesch et al. 2007; Arlt et al. 2007; Kitchatinov & Rudiger 
2007), is the presence of that large scale poloidal field which is 
allowed to evolve freely through advection by the meridional 
circulation, Ohmic diffusion, and eventually interaction with 
the instability-generated field. 

Our simulations, together with the equations of the problem 
and the resolution methods, are described in full detail in BZ06. 
We used the global ASH code (Clune et al. 1999; Brun et al. 
2004) to solve the relevant anelastic MHD equations (Eqs. 1-5 
in BZ06) in a spherical shell representing the upper part of the 
solar radiation zone (0.35 < r/R Q < 0.70), using a resolution of 
N r xN g xN^ - 193 x 128 x256. For numerical reasons, we had 
to increase substantially the diffusivities of heat, magnetic field 
and momentum, as shown in Table [3] while respecting their 
hierarchy in the solar conditions. The characteristic evolution 
times are shortened accordingly, but contrary to BZ06 we made 
no effort here to rescale them by the Eddington-Sweet time in 
order to facilitate the comparison with the real Sun. In addition 
to the case A discussed in BZ06, we performed an additional 
series of simulations with a lower Ohmic diffusivity (by a factor 
of 10, case B), in order to reach a higher magnetic Reynolds 
number. 

The temporal evolution of the magnetic fields and of the 
MHD instabilities is best followed in Fig. [3] where we display 




1x10 4 2x10 4 3x10 4 4X10 4 5x10 4 6x10 4 
Time [days] 

Fig. 3. Time evolution of the energies of the mean poloidal 
(PME), mean toroidal (TME) and non-axisymmetric (FME) 
components of the magnetic field. Cases A and B refer respec- 
tively to higher and lower magnetic diffusivity (cf. Table 0. 
Note the steady decline of the poloidal field, which is not af- 
fected by the irruption of the to = 1 Pitts & Tayler instability 
(at t^ 8,000 days in case A and » 20,000 days in case B). 



the energies of the poloidal, toroidal and non-axisymmetric 
components of the field. Initially a purely poloidal field of 
about 1 kG (when measured at the base of the computa- 
tional domain) is buried in the radiation zone; it is unstable to 
non-axisymmetric perturbations of high azimuthal wavenum- 
ber (to m 40), as shown in BZ06 (cf. their Fig. 7), which is in 
agreement with the predictions of Markey and Tayler (1973). 
The poloidal field diffuses outward at a rate proportional to 
the Ohmic diffusivity and at some point (around 8,000 days 
in case A or 20,000 days in case B) it meets the differential 
rotation that has spread into the radiation zone, due to thermal 
diffusion. Their interaction induces there a toroidal field, whose 
strength is comparable to that of the poloidal field in case A; it 
becomes even larger in case B, where it keeps growing when 
we stop the simulation, at 60,000 days (corresponding to 1.3 
Gyr when rescaled to the Ohmic diffusivity of the Sun). This 
toroidal field produces a strong non-axisymmetric MHD insta- 
bility, with a dominant azimuthal wavenumber to = 1, which 
is clearly the signature of the Pitts & Tayler instability: it is 
illustrated in Fig. |4] In case A, this instability-generated field 
saturates at an energy comparable to that of the mean poloidal 
field, whereas in case B it is still increasing when we stop the 
simulation, much like the toroidal field. 

This can be seen also in Fig. [5] where we display a merid- 
ional section of the mean and fluctuating components, from the 
case B simulation at t — 50, 000 days, where the energy of the 
fluctuating field matches that of the poloidal field. By subtract- 
ing the data of opposite azimuthal angles, <p and <p+w, we isolate 
the modes of odd to, and we choose <p such as to optimize the 
contribution of the to = 1 component. The mean field has 
a quadrupolar shape and is strongest close to the convection 
zone, whereas the fluctuating field b v is present everywhere, 
at medium level. At high latitude, the fluctuations are broader 
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Fig. 4. Azimuthal component of the instability-generated mag- 
netic field, at the peak of the instability (« 10,000 days, case 
A), measured at r/R Q = 0.7, 9 = 75°. Note that this non- 
axisymmetric field is dominated by the wavenumber m = 1, 
which is the signature of the Pitts & Tayler instability. 

in the horizontal than in the vertical direction; translated into 
wavenumbers, this anisotropy amounts to n/l * 5. Measured 
in the vertical direction, the characteristic scales are somewhat 
larger for O.O7/? ) than for (* O.O5fi ). 

7. Confronting the numerical simulations with the 
analytical model 

How do our 3-D simulations compare with the simplified ana- 
lytical model described above in §2-3? The Alfven frequency 
cda ~ 3 10~ 7 s _1 that characterizes the toroidal field at the peak 
of the instability (=s 10 4 days in case A) is significantly below 
the critical value for the Pitts & Tayler instability which is de- 
duced from ®: oJf = 2 10~ 6 yflt/IZ s -1 , where l m = In /I is 
the radial wavelength. 

Another property of the analytical model is that the vertical 
wavenumber largely exceeds the horizontal wavenumber, par- 
ticularly when the stratification is dominated by the yU-gradient, 
since (l/n) 2 (N fI /a> A ) 2 = A M = 0(1) (Eq.[l4]l. In our model, the 
contrast n/l is less severe, because there is no composition gra- 
dient: the instability threshold is then given by ( fT~8b , which we 
rewrite here as 

I 2 N 2 

s^-^-=sA t <3. (24) 

n W A 

With the parameters that characterize our numerical simula- 
tions, this yields n/l > 50, which is still not compatible with 
the observed structure of our m — 1 mode, where n/l m 5. It 
is as if the role of the stratification were overestimated in the 
analytical model. 

Therefore we must conclude that the analytical model con- 
sidered by Spruit is not applicable to the more realistic prob- 
lem treated in BZ06, and to the situation arising in stars. There 
may be several reasons for this. The first is that this simplified 
model ignores the presence of the poloidal field. However its 




-7.0e+03 7.0e+03 G 



Fig. 5. Mean toroidal field B^(r, 6) (left) and meridional sec- 
tion of the fluctuating toroidal field b v (r, 9, <p) (right). Only the 
odd m have been kept in the latter, and the azimuthal angle 
<p has been chosen such as to emphasize the m = 1 compo- 
nent. Note that mean and fluctuating fields are of comparable 
strength. (Case B, t = 50, 000 days.) 

inclusion would probably tend to stabilize the configuration, 
and that would worsen the discrepancy. 

Second, the rotation is taken as uniform, whereas the shear 
of the differential rotation certainly plays a destabilizing role. 
Moreover, it is that neglected differential rotation which is the 
very cause of the instability by generating the toroidal field, a 
point that has been stressed by Maeder & Meynet (2003). 

The third reason is that the %u dependence of the insta- 
bility is poorly described by a trigonometric function, unless 
the wavelength is small compared to the distance to the axis: 
l m <sz m. Bessel functions would probably be better suited (cf. 
Tayler 1957). 

Finally the simplified model ignores the curvature effects, 
as was emphasized by Denissenkov & Pinsonneault (2006). 
The neglect of the ?zr-component of the buoyancy force, com- 
pared to that of the Lorentz force, is allowed only as long as 
N 2 (Tu/f)^ w <K to\^ m , hence in close vicinity of the rotation 
axis, for vj/r <s: (oja/N) 2 . This is a very severe restriction, 
since (oja/N) 2 ~ 10" 6 in our numerical solution (case A). 

8. Looking for dynamo action 

An important property of our numerical solutions is that the 
decline of the poloidal field is not affected by the instability- 
generated field. As can be seen in Fig. [3] this is true even once 
the (Pitts & Tayler) instability has reached its saturation level, 
where its energy is comparable with that of the mean poloidal 
field. This has two consequences. First it proves that the small- 
est resolved scales do not act on the mean poloidal field as a tur- 
bulent diffusivity: they seem to behave rather as gravito- Alfven 
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waves, with their kinetic energy balancing their potential en- 
ergy (magnetic + buoyancy), as in the linear regime described 
by Eq. ( 1221 . Hence the saturation of the instability is not due 
to the mechanism suggested by Spruit (2002), which was in- 
spired by thermal convection, namely that the magnetic eddy- 
diffusivity rj t adjusts such as to neutralize the Pitts & Tayler 
instability, i.e. that it satisfies eq. [20] In our model - and pre- 
sumably also in stellar radiation zones - the regulation is appar- 
ently achieved through the action of the Lorentz torque on the 
differential rotation, which produces just the right amount of 
toroidal field that is required to sustain the instability, through 
small departures from Ferraro's law. 

At saturation, the mean quantities are stationary on the in- 
stability time-scale Q./o)\, which translates into the following 
condition for the azimuthal component of the momentum equa- 
tion: 



4jT P ^(m 2 Q.) = B„ 
at 



V{*njB 9 ) +<b m - V(tirfv)>* 0. 



(25) 



Here B is the mean axisymmetric magnetic field with B m being 
its meridional part, and likewise for b, the non-axisymmetric 
field generated by the instability; < . . . > designates the az- 
imuthal average. Since the characteristic scales are comparable 
for the mean and fluctuating field, this equation tells us that 
the magnetic energy of the instability field must be - crudely 
speaking - of the same order as that of the axisymmetric mean 
field, as we observe in our simulations. 

In our simulations we see no regeneration of the mean 
poloidal field: the a-effect plays a negligible role, at least up to 
the magnetic Reynolds number Rm = R 2 AQ./t] ~ 10 5 (caseB), 
for Prandtl number P m = v/77 = 1. (Evaluated with the instabil- 
ity velocity v at saturation, that Reynolds number is of the same 
order: Rm= v^/i, = (vi/cr,)/^ ~ V^n/aft/i; = t^fl/i/.) 
Note that the y8-effect, i.e. the turbulence-enhanced diffusivity, 
is also absent here; hence one should not expect much mixing 
of the stellar material. We thus conclude that in our simulations 
the Pitts & Tayler instability is unable to sustain a large-scale 
mean field dynamo, in the parameter domain that we have ex- 
plored. 

There is no sign either of a small-scale fluctuation dynamo, 
though one may argue that we inhibit this type of dynamo by 
imposing our large-scale fossil field. To check this point, we 
switched off the poloidal field at the latest stage of our low-77 
simulation (case B). As shown in Fig. [6] the toroidal field de- 
creases then rapidly, because it is no longer produced by the 
f2-mechanism, and the instability-generated field accompanies 
its decline. Thus the fluctuating field does not maintain itself, 
although the magnetic Reynolds number, Rm = R 2 AQ./j]~ 10 5 
(case B), should amply fulfill the necessary condition for a tur- 
bulent dynamo: at magnetic Prandtl number of order unity, as 
here, the critical magnetic Reynolds number is of order 100, 
according to Ponty et al. (2006). 

9. Why do our results differ from those of 
Braithwaite? 

To our knowledge, the only simulation so far that claims to 
support a dynamo operating in stellar radiation zones is that by 
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Fig. 6. Evolution of the magnetic energies after suppressing the 
mean poloidal field at t — 60, 000 days (case B). The mean 
electromotive force due to the instability-generated field pro- 
duces some amount of poloidal energy (PME, dashed line), but 
that field is too weak to prevent through O-efect the Ohmic dis- 
sipation of the toroidal field (TME, dash-dotted line), which 
causes the decline of the instability-generated field (FME, solid 
line). There is no dynamo action, in spite of the high magnetic 
Reynolds number Rm=10 5 . 



Braithwaite (2006): he showed that a sufficiently strong differ- 
ential rotation can amplify a seed field to a level where it seems 
to be maintained, while undergoing cyclic reversals. According 
to him, his results confirm the analytical expectations of the 
role of the Pitts & Tayler instability, but to us it is not clear 
whether that specific instability plays any role in his simula- 
tion: for instance the author does not mention the m = 1 sig- 
nature of the instability-generated field. Instead, he may have 
triggered a fluctuation dynamo. 

Why do we reach different conclusions about the existence 
of such a dynamo? Our set-ups differ somewhat, even when we 
suppress the large-scale poloidal field: in Braithwaite's case, 
differential rotation is enforced by a body force with strong re- 
laxation, whereas in ours it spreads from the boundary of the 
computational domain, which is more realistic. 

We differ also in the boundary conditions applied to the 
magnetic field, which are known to play a sensitive role in 
numerical dynamo. We connect our internal field to a poten- 
tial field outside, as if the convection zone were a perfect con- 
ductor, whereas Braithwaite imposes the field to be normal to 
the boundary. Using a geometry similar to that of Braithwaite, 
Gellert et al. (2007) do not find dynamo action either, and they 
loose the Pitts &Tayler instability when they switch off the ex- 
terior field, much as we find when we switch off the poloidal 
field (Fig©. 

But the main difference perhaps resides in the way the 
equations are solved. Our code uses (enhanced) physical dif- 
fusivities; it is of pseudo-spectral type with a resolution of 
128x256x192, and this method is known to have exponential 
convergence and machine accuracy in evaluating derivatives. 
This allows us to reach a magnetic Reynolds number of 10 5 , 
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and when we fail to observe dynamo action this is certainly not 
due to a insufficient resolution. Braithwaite uses instead a 6th 
order finite difference scheme, with a resolution of 64x64x33; 
the numerical diffusion is tuned to ensure stability for the cho- 
sen resolution, but it is not straightforward to infer from it the 
magnetic Reynolds number that characterizes the simulation 
(Braithwaite, private communication). 

Further examining his numerical code (Nordlund & 
Galsgaard 1995), one is led to guess that the background 
viscosity characterizing the numerical dissipation is v = 
0.02 c s Ax, with c s being the sound speed and Ax the grid spac- 
ing. This translates into a viscous time, which is also the Ohmic 
decay time, T n — r v = L 2 /v = 1600 r s , with r s = L/c s . It ap- 
pears that none of the simulations reported in his paper has 
been run for more than 3000 t s ; so one may wonder whether 
they have been carried long enough to go beyond the transient 
phase. 

10. Conclusion 

We have re-examined the non-axisymmetric instabilities affect- 
ing a toroidal magnetic field in a rotating star, which have 
first been described by Pitts and Tayler (1985) in the ideal, 
non-dissipative limit. The problem was generalized by Spruit 
(1999) to include the diffusion of heat and of magnetic field. 
We have extended his analytic treatment to the case where the 
medium is stratified both in entropy and in chemical composi- 
tion (Eq.r20b. Our exact solutions fully validate his approximate 
results. 

Then we have compared these analytical results with nu- 
merical solutions built with the 3-dimensional ASH code; in 
our model the toroidal field is produced by shearing a fossil 
poloidal field through the inward propagating differential rota- 
tion imposed by the convection zone. Our numerical solutions 
clearly display the Pitts & Tayler instability with its dominant 
m — 1 mode, but they do not conform to the quantitative pre- 
dictions of the analytical model. In our simulations the instabil- 
ity occurs well below the threshold predicted by the analytical 
model, and it is much less sensitive to the stratification. These 
discrepancies are probably due to the approximations made to 
simplify that analytical model, such as neglecting the poloidal 
field, the differential rotation and the radial component of the 
buoyancy force. 

It also appears that the saturation of the instability can- 
not be ascribed to a turbulent diffusivity fulfilling the critical 
conditions, as in Spruit's model: it occurs when the energy of 
the instability-generated field reaches approximately that of the 
mean fields. The mean poloidal field steadily declines due to 
Ohmic dissipation, while it is wound up by the differential rota- 
tion to produce the toroidal field. Contrary to Spruit's expecta- 
tion, which is based on questionable grounds as we have shown 
in $8] we detect here no sign of a dynamo that could regenerate 
the mean fields; the small scale motions do not act either as an 
eddy diffusivity on the mean poloidal field. Unlike the turbulent 
motions present in a convection zone, the instability-generated 
motions produce here no a and no /J-effect. Neither do we 
observe a fluctuation dynamo, in spite of the relatively high 
magnetic Reynolds number, in contrast with the findings of 



Braithwaite (2006), who however considers a somewhat sim- 
pler model. 

But the Pitts & Tayler instability persists as long as the 
toroidal field remains of sufficient strength, i.e. a few Gauss 
in the conditions prevailing below the solar convection zone, 
which puts a similar requirement on the poloidal field. We have 
shown in Brun & Zahn (2006) that such a poloidal field does 
not exist in the Sun, because it would imprint on the radiative 
interior the differential rotation of the convection zone, and that 
is ruled out by the helioseismic diagnostic. The Pitts & Tayler 
instability could well occur in other stars hosting a large-scale 
toroidal field, but we doubt that it may cause there any signif- 
icant transport of matter and angular momentum, since in our 
simulations the motions associated with the instability behave 
rather as Alfven waves than as turbulence. To settle that issue, 
observational tests will play an irreplaceable role. 
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Appendix A: Dispersion relation and energy 
equation 

Let us recall that all perturbations (displacement magnetic 
field b, pressure P', temperature T', molecular weight //) are 
expanded in Fourier modes 



exp iilm + imp + nz — cri), 



(A.l) 



in the vicinity of the rotation axis (the z-axis). 

We begin by deriving the buoyancy force, which is weak- 
ened through radiative and atomic diffusion. We split the buoy- 
ancy frequency in two parts, the first due to the thermal strati- 
fication and the second to the the composition gradient: 

N 2 = N? + Nl = -^(V ad - V) + J- 1 ^t) (A.l) 
' » Hp Hp \dlnP) 

with the usual notations, and taking for simplicity the perfect 
gas equation of state. 

The linearized heat equation may be written as 

ii 



V Nf 
If — -\ (—lO~l; z ) 



(A.3) 



where k is the thermal diffusivity, and s 2 = I 2 + m 2 jxn 2 + n 2 . 
(We simplify the Laplacians by assuming that the diffusivities 
do not vary much over a meridional wavelength.) Likewise, the 
advection/diffusion equation for the molecular weight pertur- 
bation takes the form 



-As 2 ^, 



(A.4) 

P- g P 
with A being the molecular diffusivity. Thus the buoyancy force 
is given by 



P 

1 

' P 



r 
Y 



Nf 



m 



1 + iks 1 I a- 1 + i As 2 jcr 



ft. (A.5) 



Turning next to the Lorentz force, we perturb the toroidal 
magnetic field B, = B(m) by the displacement and draw 
the field perturbation b from the induction equation 



— = Vx(VxB)-Vx (77V x B), 

ot 



to obtain 

- icrb = V x (-/cr£ x B t ) - rjs 



2 b. 



(A.6) 



(A.7) 



First ignoring Ohmic diffusion, we get b = im (B / m) It is 
then straightforward to calculate the perturbation of the Lorentz 
force per unit volume 



— [(V x b) x B + (V x B) x b] ; 

Anp 



(A.8) 



introducing the Alfven frequency u) 2 A - (B 2 lm 2 )IAnp and as- 
suming that B oc %n p \ 



L m = -m co A % m + (Im) mu> A ^ ip - im (p + 1) u> A ^ 

Lip = im(p + \)u? A % m 

L z — (nvj) mti) A ^,p — m 2 a>A 



(A.9) 



When keeping Ohmic diffusion b = im(B/m)g/(l + irjs 2 /cr), 
and hence the Lorentz force will also be divided by ( 1 + /77 s 2 /cr). 

It remains to implement the expressions derived above for 
the buoyancy and Lorentz forces in the equations of motion: 



i P ' 

il — + 

P 



2 2 

m a) . 



1 + irjs 2 /cr 



(lm)mio 2 A (p + l)mu> A 



+ 1 



1 + irjs 2 l<r 1 + irjs 2 I<t 



. m P' 

■uj p 

P' 
in 



(p + \)moj 2 A 
1 + irjs 2 l<r 

'I 



+ 2;Qo- 



o 2 ^ 



(nvj)m<ji) 2 



p 1 + irjs 2 jcr 



-a 2 + 



2 2 
m lo a 



Nf 



& = 0. 

N 2 , 



1 + ir\s 2 jcr 1 + ixs 2 Iyer 1 + i As 2 /cr 



(A. 10) 



and to complete them with the continuity equation, in the 
Boussinesq approximation: 



l£ m + —ft, + n& = 0. 



(A.ll) 



This yields a fourth-order system whose determinant is the dis- 
persion relation we are looking for. 

From here on, we shall consider only the most realistic case 
where p — 1 . To first approximation we may neglect the molec- 
ular diffusivity compared to the Ohmic diffusivity, and a fortiori 
to the thermal diffusivity. We further assume that I 2 «: n 2 , scale 
all frequencies and damping rates by the Alfven frequency: 



cr 

cr = = a + ip, 



Q = — , 



Kfl 

ya)A 



rjn 



and introduce 

I 2 N 2 
A - — -L 

' ~ 2 2 ' 



7 2 N 2 

A - —-± 

f ~ 2 2 ' 



^, (A12) 



(A.13) 



The dispersion relation may then be cast into 

cf 6 - c? 4 [4Q 2 + A, + Af, + 2 + 2hk + h 2 ] - d 3 8Q/w 

+ a 2 [A,+ A^-3 + 2hk(AQ}+ A^ + 1) + /i 2 (4Q 2 + A, + A^)] 

+ crhk8£l/m - hkA^ 

+ i <x 5 [2h + k] 

- ia 3 [k(4Q 2 + A tl + 2) + 2h (4Q 2 +A,+A I1 + Y) + h 2 k] 

- ia 2 8(k + h)Q./m 

+ ia [k (A/j - 3) + h (A, + A^) + h 2 k(4-h 2 + A^)] = (A. 14) 

From now on we no longer deal with the gravito-inertial 
modes, which allows us to discard the terms in <x 6 and cr 5 . The 
slow modes then obey a 4th order equation which can be further 
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reduced to a third order equation, by taking into account that 



4Q 2 » (sA„Au), s « 1 and cr 2 <k (l,k 2 ). This leads to 



- iS£(£ 2 +4£ + 3-A*-S 2 ) = 0, (A.15) 

where £ = 2mQxr = 2m£l <t/u> a , 5 = 2Q/j = 2Qrin 2 /a> A and 
A* = eA, + A/j. 

This dispersion relation may be solved numerically to ob- 
tain the growth-rate 07 of the unstable modes. In the two limit 
cases, this equation has solutions for S — > at given A^ or sA,, 
which can be easily derived. When A, - 0, the growth-rate is 
given by 



, 4 JTTJ^ - 2 - A„ 

o-i = T]n 2 —!- — ■ (A, = 0, S -> 0); (A.16) 

2 + 2A,, - 4 yi + A p 

it is positive for A^ < 6 +4 a/3. In the other special case A M = 0, 
the growth-rate is positive between sA, = and 3: 

cr 7 = rjn 2 3 ~f A ' (A„ = 0, S -> 0). (A.17) 
eA, 

These expressions are not valid for A^ — > 3, or respectively 
A, — > 0, eA ( — > 3, where one has to keep higher order terms. 

We finally turn back to the non-dissipative case, for which 
we derive the dynamical energy equation. Introducing the per- 
turbation velocity v = dg/dt, and considering only the real part 
of all variables, we rewrite ( IA. lOt for the stable case p = 1 as 

— -2Qv^ + - — 
at p dm 

+m 2 aj 2 A ft, - m(lm) cj 2 a ft, + 2 o? A ^ = (A. 1 8) 

^ +2Qv . + J-f:- 2.^=0 (A.19) 
at pm o<p o<p 

dv 1 BP' 

— - +- — + (N 2 + m 2 o? A ) ft - m{nm) u\ ft = 0. (A.20) 
at p oz v ' 

We multiply these equations respectively by v m , v v , v z and add 
them up, making use of the continuity equation ( IA.1U . to ob- 
tain 

I A \ v 2 + v 2 + v 2] + I m 2 2 d_ r 2 + £ + gl 1^2 ^2 

2 dt*- v z i 2 A dt^ m tz J 2 ar* 



2^ 



<9ft^<9ft _ dft dft, 
df d(p dt dip 



(A.21) 



For these stable oscillatory modes, the r.h.s. vanishes since 
d/dt = icr and d/d<p = im; therefore the sum of kinetic and 
potential (magnetic + buoyancy) energies is constant, as is the 
rule in non-dissipative dynamical systems: 

\v 2 + l -m 2 oj 2 A f + l -N 2 f z = est. (A.22) 



